Universal thermal and electrical transport near the 
superconductor-metal quantum phase transition in nanowires 
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We describe the thermal (k) and electrical (a) conductivities of quasi-one dimensional wires, across 
a quantum phase transition from a superconductor to a metal induced by pairbreaking perturbations. 
Fluctuation corrections to BCS theory motivate a field theory for quantum criticality. We describe 
deviations in the Wiedemann-Franz ratio k/oT (where T is the temperature) from the Lorenz 
number {tv'^ /3){kB /e)^ , which can act as sensitive tests of the theory. We also describe the crossovers 
out of the quantum critical region into the metallic and superconducting phases. 



The Wiedemann-Franz law relates the low tempera- 
ture (T) limit of the ratio W = K,/{aT) of the thermal (k) 
and electrical (a) conductivities of metals to the univer- 
sal Lorenz number Lq = (7r^/3)(fcB/e)^. This remarkable 
relationship is independent of the strength of the inter- 
actions between the electrons, relates macroscopic trans- 
port properties to fundamental constants of nature, and 
depends only upon the Fermi statistics and charge of the 
elementary quasiparticle excitations of the metal. It has 
been experimentally verified to high precision in a wide 
range of metals and realizes a sensitive macroscopic 
test of the quantum statistics of the charge carriers. 

It is interesting to note the value of the Wiedemann- 
Franz ratio in some other important strongly interact- 
ing quantum systems. In superconductors, which have 
low energy bosonic quasiparticle excitations, a is infinite 
for a range of T > 0, while k is finite in the presence 
of impurities Q, and so W — 0. At quantum phase 
transitions described by relativistic field theories, such as 
the superfluid-insulator transition in the Bose Hubbard 
model, the low energy excitations are strongly coupled 
and quasiparticles are not well defined; in such theories 
the conservation of the relativistic stress-energy tensor 
implies that k is infinite, and so = oo Q. Li and 
Orignac [3] computed W in disordered Luttinger liquids, 
and found deviations from Lq, and found a non-zero uni- 
versal value for W at the metal-insulator transition for 
spinless fermions. 

The present paper will focus on the quantum phase 
transition between a superconductor and a metal (a 
SMT) . We will consider quasi-one dimensional nanowires 
with a large number of transverse channels (so that the 
electronic localization length is much larger than the 
mean free path (iY) which can model numerous recent 
experiments 0, 0,0, i, S, [H [HI, [H, El. We will de- 
scribe universal deviations in the value of W from Lq, 
which can serve as sensitive tests of the theory in future 
experiments. 

The mean-field theory for the SMT goes back to the 
early work [l4| of Abrikosov and Gorkov (AG): in one 
of the earliest discussions of a quantum phase transition, 



they showed that a large enough concentration of mag- 
netic impurities could induce a SMT at T = 0. It has 
since been shown that such a theory applies in a large 
variety of situations with 'pair-breaking' perturbations: 
anisotropic superconductors with non-magnetic impuri- 
ties [l5| . lower-dimensional superconductors with mag- 
netic fields oriented in a direction parallel to the Cooper 
pair motion [l6| , and s-wave superconductors with inho- 
mogeneity in the strength of the attractive BCS inter- 
action 17|. Indeed, it is expected that pair-breaking is 
present in any experimentally realizable SMT at T 0: 
in the nanowire experiments, explicit evidence for pair- 
breaking magnetic moments on the wire surface was pre- 
sented recently by Rogachev et al. 



Fluctuations about the AG theory have been con- 
sidered [16i, IS, 19] in the metallic state, and lead to 
the well-known Aslamazov-Larkin (AL), Maki-Thomson 
(MT) and Density of States (DoS) corrections to the con- 
ductivity. At the SMT, field-theoretic analyses ^ 
show that the AG theory, along with the AL, MT and 
DoS corrections, is inadequate in spatial dimension d < 2, 
and additional self-interactions among Cooper pairs have 
to be included. Here, d defines the dimensionality of the 
Cooper pair motion, while the metallic fermionic quasi- 
particles retain a three-dimensional character; therefore, 
the confining dimension, R, is larger than the inverse 
Fermi wavevector, but smaller than a superconducting 
coherence length or Cooper pair size, ^. The behavior 
of W has been considered in this field-theoretic frame- 



work 21|, and it was found that there were logarithmic 
corrections to the Lorenz number in d = 2. Here we 
will examine the d = 1 case in some detail: the tran- 
sition is described by a strongly-coupled field theory of 
bosonic Cooper pairs, overdamped by their coupling to 
the fermionic quasiparticles. Remarkably, all important 
couplings between the bosons and the fermions scale to 
universal values, and consequently the Wiedemann-Franz 
ratio of this theory also approaches a universal constant 
which we compute in a 1 /N expansion (the physical case 
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FIG. 1: (Color online) Crossover phase diagram of the 
superconductor-metal transition in a quasi-one dimensional 
superconductor. The "Metal" is described by the perturba- 
tive theory of Ref. [l6l . The "Quantum critical" region is de- 
scribed by S and realizes our result for W in Eq. ([TJ. The 
Mooij-Schon mode is present everywhere, but couples strongly 
to superconducting fluctuations only in the "Fluctuating su- 
perconductor" regime, where it is described by Eq. Q; note 
that S does not apply here. The dashed lines are crossover 
boundaries which, by Eq. Q occur at T ~ |a — ac\^^ . 



systems. Our theory includes amplitude and phase fluc- 
tuations on an equal footing, along with strong damping 
from the fermionic modes, and describes the transition 
into a metallic phase at T = 0: this is the case even 
though physics of the Mooij-Schon mode is present in 
the "fluctuating superconductor" regime in Fig. [TJ 

It is useful to place our results in the context of re- 
cent microscopic computations in BCS theory [l6| on the 
metallic side, with the pairbreaking parameter a larger 
than critical ac of the SMT. These results were obtain in 
the dirty limit (£ <C '^), but we expect the same theory 
(the action S below) to apply in the quantum critical 
regime in the both the clean and dirty limits (although 
there are distinctions in the "fluctuating superconduc- 
tor" regime of Fig.[T|). For the conductivity, these results 
are lldl 
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Our present computation of W ignores the influence of 
disorder on quantum criticality, and this may require the 
clean limit ^ < £ [111. 

We also discuss the nature of the crossovers from this 
universal quantum critical physics to previously studied 
regimes at low T about the superconducting and metallic 
phases: these are summarized in Fig. [TJ On the metallic 
side, there is a crossover to a low T regime described by 
the theory of AL-l-MT-l-DoS corrections in d = 1. On 
the superconducting side, there is a regime of intermedi- 
ate temperatures where the classical phase slip theory of 
Langer, Ambegaokar, McCumber, and Halperin (LAMH) 
applies 23|, |2J| , and eventually another crossover at still 
lower temperatures to a phase fluctuating regime whose 
description requires a non-linear cr-model of fermion pai r 
fluctuations coupled to the superconducting order [25| . 
Here the phase fluctuations are essentially equivalent to 
a plasma charge oscillation, which in d = 1 is the Mooij- 
Schon mode [ig] . A number of works [i^, [H, 0, [sO] have 
examined the destruction of superconductivity due to 
quantum phase slips in such a phase fluctuating regime: 
we maintain that the phase with phase slip prolifera- 
tion in such models is an insulator at T = 0, and so 
such theories describe a superconductor-insulator quan- 
tum transition, and may be appropriate in inhomogenous 



where ctq is a background metallic conductivity, c is a 
non-universal constant, D is the diffusion constant in the 
metal, and the remaining corrections from pairing fluc- 
tuations have been written in the form of a power of T 
times a factor within the square brackets which depends 
only upon the ratio h{a — ac)/kBT. This way of writing 
the results allows us to deduce the importance of the fluc- 
tuations corrections, in the rcnormalization group sense, 
to the SMT. The flrst square bracket represents the usual 
AL correction, and has a prefactor of a negative power 
of T, and so is a relevant perturbation; this is so even 
though this correction vanishes as T ^ 0. The second 
square bracket arises from the additional AL, MT and 
DoS corrections: the prefactor has no divergence as a 
power of T, and so this correction is formally irrelevant 
at the SMT. Note, however, the complete second term 
has a finite limit as T — > 0, and so becomes larger than 
the formally relevant AL term at sufficiently low T in 
the metal. We therefore identify the second term as dan- 
gerously irrelevant in critical phenomena parlance: i.e. 
important for the properties of the low T metallic re- 
gion, but can be safely neglected in the shaded quantum 
critical region of Fig. [TJ 

Armed with the above insights, we focus on the fiuc- 
tuations associated with the usual AL correction. These 
have [3] a Cooper pair propagator {Dq'^ + \lu\ + a)~^ 
at wavevector q and imaginary frequency cu in the metal 
in both the clean and dirty limits. This motivates the 
quantum critical theory of Ref. j2Q for a field '^{x, r) rep- 
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resenting the local Cooper pair creation operator: 



Eq. Q is 



S 



dx 



dr 



(3) 



This theory will apply to quasi-one dimensional wires for 
R < ifiD/kBTY''^. In the dirty limit, ^ < 

we have 

D = D, but the value of D is different in the clean case. 
All couplings in S are random functions of position; in 
particular, randomness in a is expected to be relevant at 
the quantum critical point. We neglect this randomness 
in our quantitative results in the quantum critical region, 
and so they only apply above a T ~ Tjis which can be 
made arbitrarily small in the clean limit. 

From S, we obtain[20] the singular contribution to the 
conductivity in the vicinity of the quantum critical region 
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where z is the dynamic critical exponent, is a correla- 
tion length exponent, and <i>a- is a universal scaling func- 
tion. In a Gaussian approximation, the Kubo formula 
yields z = 2, y = 1/2, and $<,(y) = [tt / {I2^))y~^'^ , 
and so this result is in precise correspondence with the 
first term in Eq. ([2]). In the limit T <C (a— Oc)^'', we have 
already seen that this term is subdominant to the danger- 
ously irrelevant second term in Eq. However, moving 
into the quantum critical region where T ^ {a — acY'^ ^ 
the contribution from Eq. (|4]) dominates all other terms, 
and we have a ^ T^^/^<I>o-(0). The microscopic analysis 
of Ref. ^ obtained a ^ T^^, which is valid only at T 
large enough where u can be neglected. Going beyond 
the Gaussian theory, the values of z and v have been 
determined in a c? = 2 — e expansion l2Ct ISl | , and also 
in quantum Monte Carlo simulations |32| with excellent 
agreement. Here, we have obtained these exponents in a 
theory with N complex fields ^ directly in d = 1 , and 
obtained to order 1/iV 



0.131 

N 



V = 1- 



0.389 



(5) 



to be compared with the Monte Carlo estimates of z = 
1.97 and v = 0.689 [11]. The value of $^(0) has a non- 
universal cutoff dependence associated with anomalous 
dimension 2 — z. Note that the results above for a in 
the metallic and quantum-critical regimes imply a non- 
monotonic T dependence for a > ac, possibly consistent 
with the observations of Ref. 
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Similar reasoning can be applied to the thermal con- 
ductivity K, which can be computed from S using a sep- 
arate Kubo formula [s^. The scaling form analogous to 
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with another universal function. We have verified 
that the Gaussian prediction from S again agrees with 
the perturbative AL contribution of the microscopic the- 
ory at low T. Our main result for the Wiedemann-Franz 
ratio in Eq. ^ follows from W = (fcB/e)2$K(0)/$^(0), 
as the nonuniversal prefactor cancels out in ratio of these 
scaling functions; the 1/N expansion was carried out by 
generalizing the methods of Ref. 's?. 

We now turn to an important conceptual issue: the 
role of charge conservation and associated normal modes. 
From hydrodynamic arguments we know that a one- 
dimensional metal or superconductor should support a 
gapless plasmon, or a Mooij-Schon normal mode [l^. 
In d = 1, this mode is gapless and disperses as ^ 
q\n^^\l/{qR)). In our theory for the quantum critical 
region, the Cooper pair field carries charge 2e but only 
exhibits diffusive dynamics with to ~ bq^, and there is no 
Mooij-Schon mode in the dynamics of the action S. The 
answer to this puzzle is contained in arguments made 
in Refs. 35 and [s^ on the role of conservation laws in 
the critical fluctuations of quantum transitions in metal- 
lic systems for which the order parameter is overdamped 
(as is the case here). These early works considered the 
onset of spin-density wave order in a metal; in the quan- 
tum critical region, the spin excitations consisted of dif- 
fusive paramagnons whose dynamics did not conserve to- 
tal spin. However, loffe and Millis [1^ argued that the 
Ward identities associated with spin conservation only 
imposed significant constraints on the effective action at 
uj ^ q, and played little role in the w ~ regime im- 
portant for the critical fluctuations. Essentially the same 
argument can be applied here: the Mooij-Schon mode is 
present only at relatively high frequencies w ~ and the 
critical theory S describes the overdamped Cooper pair 
modes in the distinct region of phase space with ~ g^. 
The Mooij-Schon fluctuations lead to oscillations in the 
local electrochemical potential, but these remain essen- 
tially decoupled from the critical modes described by S 



3a | (see however, Eq. ^ below). It must be noted 
that the action S is not valid for lo ^ q and a complete 
description in terms of the underlying fermions is neces- 
sary to obtain the proper dynamics, which will contain 
the Mooij-Schon mode, as required. 

Further insight into this issue is gained by lowering the 
temperature from the quantum critical regime into the 
"fluctuating superconductor" regime of Fig.[T]for a < ac- 
When ksT < [ac — a) the action S does not apply for 
the smallest wavevectors and frequencies. The reasons 
for this are again analogous to arguments made for the 
spin-density-wave ordering transition in metals, as dis- 
cussed in Ref. |37|. For the latter case, it was argued that 
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with the emergence of long-range spin density wave order, 
the low energy fermionic particle-hole excitations at the 
ordering wavevector were gapped out, and so the diffu- 
sive paramagnon action applied only for energies larger 
than this gap. At energies smaller than the gap, spin- 
waves with dispersion to ^ q emerge, as a consequence 
of Ward identities associated with spin conservation. In 
the superconducting case of interest here, there is no true 
long-range superconducting order at any T > 0, but the 
order is disrupted primarily by 'renormalized classical' 
thermal fluctuations of the phase, of the complex 
field. We can now assume that there is a local pair- 
ing amplitude in the fermion spectrum, analogous to the 
spin-density wave order. Charge conservation plays an 
important role in the effective action of 4> fluctuations, 
and the 'spin-waves' in this case are, of course, just the 
Mooij-Schon excitations. The low energy effective action 
for (j) cannot be obtained from 5; rather, we have to inte- 
grate the fcrmions out in the presence of a local pairing, 
as outlined above ^7|, and we expect an action 



27r J 27r41n(l/((7i?)) 
+ I dxdr Ki{dr(l) - 2eArf + K2{d^<j)f] (7) 



For the spin density wave case [36l | it was found that 
Ki ^ (ttc — a)^/^ and K2 ~ (etc — a), and we expect sim- 
ilar behavior here. We have explicitly included the action 
of the scalar potential iAr which mediates the Coulomb 
interaction in one dimension. The normal modes of 5^ 
are the Mooij-Schon oscillations which arc now identified 
with the fluctuations of the superconducting order. This 
should be contrasted from the situation in the shaded 
quantum critical region of Fig. [1] where the Mooij-Schon 
oscillations were decoupled from the critical modes de- 
scribed by S. 

Actually, there is a weak coupling between the criti- 
cal modes of ^I* in 5 and the Mooij-Schon mode due to 
particle-hole asymmetry (there is no analog of this phe- 
nomenon for the spin density wave case). In general, the 
critical 5* action contains the perturbation [20| 



Si 



dxdr 



7** ( ^ - 2eiAr ] * 
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where 7 is proportional to the energy derivative of the 
density of states, and hence small. Indeed, the same ratio 
of the pairing to Fermi energy, which justified the present 
quasi-one dimensional treatment of Cooper pairs but not 
electrons, also causes 7 to be small. We examined the 
renormalization group fate of the perturbation Si at the 
fixed point of S: we computed the two-loop fiow of 7 in 
the d — {2 ^ e) expansion and found 



So 7 is relevant, but the scaling dimension is extremely 
small. Along with small bare value of 7, such particle- 
hole asymmetric effects can justifiably be ignored in ex- 
perimental applications. 

The main experimentally testable results of this pa- 
per are: the crossover phase diagram in Fig. [l] the 
Wiedemann-Franz ratio in Eq. ([T]) which applies in the 
shaded region, and the non-monotonic T dependence of 
a for a > ac- 
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